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ABSTRACT 

With recent advances in microcomputer technology, on board data reduction for 
orbit determination is now feasible. On-board data processing would be advan- 
tageous for a class of satellite applications which involve a one way Doppler posi- 
tioning system. Examples are the Nimbus-RAMS system, satellite aided search 
and rescue, and coastal surveillance in support of the newly imposed 200 mile 
fishing limit. Orbit determination requirements for these systems are typically 
in the region of 1 km. 

This paper presents a proposal for on-board orbit determination which would rely 
entirely on reference beacons located within continental United States. The bea- 
cons would emit a one way Doppler signal in which is encoded coordinates of the 
beacon’s location. An on-board computer would process the signal and update 
the satellite ephemeris. When a Doppler signal external to the system is received, 
the problem is inverted in the sense that the satellite ephemeris is considered as 
known and the coordinates of the Doppler signal source are estimated and trans- 
mitted to the ground. The feasibility of this concept is dependent on the possi- 
bility of developing algorithms which are at once compatible with limitations of 
on-board computers and accurate enough to maintain an adequate satellite 
ephemeris. 

A conventional batch-processing filter is not ideal for this application because of 
its data storage requirements. Another problem associated with batch-processing 
of satellite tracking data is that errors in the models of satellite dynamics are such 
that in order to achieve acceptable accuracies the time span of the data must be 
limited. A recursive filter appears to be a more logical choice. But conventional 
recursive processing schemes such as the usual form of the Kalman filter, which 
do not explicitly account for dynamic modeling errors, encounter similar diffi- 
culties. In this case, the problem manifests itself through the computed covari- 
ance matrix of the estimate reducing towards zero. This in turn results in new 
measurement information being effectively ignored. Consequently, measurement 
residuals will increase, manifesting the well-known Kalman filter divergence 
phenomenon. 
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The filter used to generate the results of this paper is recursive from one pass of 
Doppler data to the next. However, each pass is processed in a conventional batch 
mode. A “state noise” covariance matrix is added to the computed covariance 
matrix of the estimate at each step in the recursion, thus preventing the eigen 
values of the matrix from approaching zero. The form of the state noise covari- 
ance matrix is chosen from a priori mathematical considerations and then modi- 
fied by a multiplicative constant obtained from an examination of residuals. The 
resulting algorithm is sufficiently simple that it can be implemented on a small 
computing machine. 

Numerical simulations were performed to test the validity of the procedure. 
Reference beacons were situated.at Portsmouth, New Hampshire and San Fran- 
cisco, California. Epoch elements for a typical TIROS-N orbit were utilized. 

Data was generated using a 1 3 x 1 3 degree reference field. A random number 
generator added white noise of standard deviation 1 m/sec. The data was reduced 
using a modified 4x4 degree field. Maximum observed position errors were less 
than 1 km. The same procedure was implemented for real data obtained from 
Nimbus reference beacons with comparable results. 
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Introduction 

With recent advances In the technology of 
low power, low weight computers, on-board data 
reduction for orbit determination is now 
feasible [l,2]. On-board data processing would 
be advantageous for a class of satellite appli- 
cations which involve a one-way Doppler posi- 
tioning system. Examples are the NIMBUS -RAMS 
system [3], satellite aided search and rescue 
[4,5], and coastal surveillance in support of 
the newly imposed 200 mile fishing limit [6], 

Orbit determination accuracy requirements for 
these systems are typically in the region of 
1 km. 

An examination of the NIMBUS-RAMS system 
provides insight into the possible usefulness 
of on-board orbit determination and data pro- 
cessing. Typically, 10 to 15 one-way Doppler 
measurements from each of up to 50 sources are 
recorded on the NIMBUS-6 and telemetered on 
command to a control station such as NASA- 
Fairbanks, Alaska. The data is then relayed 
to a central computer at the Goddard Space 
Flight Center. With the aid of ephemeris 
information obtained from independent tracking, 
the data is processed and the results trans- 
mitted to individual users. The delay time is 
several days. The need for faster response 
has been expressed by several members of the 
user community. Also, future applications of 
the satellite Doppler positioning concept 
(coastal surveillance, satellite aided search 
and rescue) will require near real time signal 
processing. 

The need for rapid signal processing could 
be satisfied if tracking data were processed on- 
board a satellite. In that case ephemeris infor- 
mation based on fresh tracking data would be 
available at all times to the satellite's 
computer* Hence a Doppler signal received by 
the satellite could be processed by the on-board 
computer and the resulting position estimate 
transmitted to the ground. Very low telemetry 
data rates are adequate to transmit position 
information. With low data rates, the use of 
omnidirectional or low-gain, wide-beam-width 
antennas becomes possible and costs are reduced 
by simplifying the antenna array and antenna 
drive assembly. Low data rates also permit the 
use of a geostationary satellite system such as 
the Tracking and Data Relay Satellite System [7] 
to relay information to inexpensive receiver 
terminals. Another significant fact is that the 


number of Doppler signal sources monitored per 
satellite pass can be significantly increased as 
a consequence of the relatively low bit alloca- 
tion required for each position determination. 

This paper presents a proposal for on-board 
orbit determination which would rely entirely on 
reference beacons located within the Continental 
United States. The beacons would emit a signal 
in which is encoded coordinates of the beacon's 
location. An on-board computer would process the 
signal and update the satellite ephemeris. When 
a Doppler signal external to the system is re- 
ceived, the problem is inverted *in the sense that 
the satellite ephemeris is considered as known and 
the coordinates of the Doppler signal source are 
estimated and transmitted to the ground. The 
feasibility of this concept is, of course, depen- 
dent on the possibility of developing algorithms 
which are at once compatible with the limitations 
of on-board computers and accurate enough to main- 
tain an adequate satellite ephemeris. This paper 
recommends a filter which is recursive on a pass- 
by-pass basis and which has a fading memory to 
account for the degrading effect of gravity field 
error. The recursive filter requires the storage 
of only one pass of data at a time, and hence is 
suitable for on-board computers with limited 
storage capacity. The mathematical details are 
described in the succeeding sections. The paper 
also provides the results for simulated data and 
real data reductions as well as suggestions for 
further analysis and development of these 
concepts. 

A Recursive Filter for Satellite State Estimation 

Let {Z(Ti)}i£$ be a set of Doppler passes 
with each pass referenced to a time Ti which may 
be taken as the time for the first observation 
in the pass. Assume that observation set 
{Z(Ti) )i^n has been processed to yield a minimum 
variance estimate, 6x(T^) , of the deviation, 
6X(Tfl) , of the satellite state at time Tn from a 
nominal orbit. Data pass Z(T$+i) is to be com- 
bined with 6 X(Tjj) in a way which yields a^suf f i- 
ciently accurate estimate, <5X(T$+i), of 6X(Tfj+j) , 
Furthermore, the resulting algorithm is to be 
compatible with small computing machines. 

Assume a linear state transition equation 

6X(T N+1 ) - <K%h, %) 5X(T n ) + T m (I) 

where $(%+]., t n) is a six-by-six state transi- 
tion matrix and where Tn+ 1 is a six dimensional 
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state noise random vector. Equation 1 is 
assumed to be referenced to a standard inertial 
coordinate set. Define 


Under assumptions A through D, one can write 
Tn+1 = (10) 


E(6X(T N ) 6xT(T N )) = P N (2) 

Before the data set, Z(Tj|+l) , is processed, 
the optimal estimate of 6 X( Tjb- 1) is obtained 
from Equation 1. 

dX-CTN+i) - (f>(T N4 -l,T N ) <5X(T N ) (3) 

Hence 

e(6X-(T n+ i) 6X-(T N+ 1)) = P N+1>0 = 

< t>( T lH-l» T N) PN^ T ( T N+1. T N) = ^N+l 

where 

e ^ t N+1 tT N+1^ = AP N+1 

Define a linearized observation equation as 

6Z(X N+ 1) = A N+1 6X(T N+1 ) (6) 

where 6 Z(Tn+i) represents a vector of deviations 
of the noiseless or correct representations of 
the data from nominal values. Since the data are 
corrupted by noise, we have 

e((6z(X n+1 )-6z(T n+1 ))(6z(T n+1 )-6z(T n+1 ))^= q (7) 

Given Equations 1^ through 7, a minimum 
variance estimate of 6 X(Tn+ 1) is 

fiXCTN+i) = (A.TQ" 1 A+Pn+! , q) " 1 ( A t Q" 1 6 Z (T n+1 ) 

(8) 

+ P S-l,0 6r(T N+l)) 


and 

E((6X(T N+1 )-6X(T N+1 )) (6X(T n+1 )-6X(T n+1 ))t) = 

(9) 

P N+1 = (A T Q -1 A+Pn| 1>0 ) _1 

It remains to specify how the covariance 
matrix, APjj+x °f T N+1> is obtained. To generate 
the results of this paper, the following assump- 
tions are imposed: 

A. Components of are uncorrelated 

in the instantaneous along-track, cross- 
track, and radial coordinate set. 

B. Along-track, cross- track, and radial 
components of tn for both position and 
velocity are in the ratio of 10:5:2. 

C. The along-track position component 
of Tfl-i-i is a linear function of Tji+j-Tfl, 

D. The along-track velocity component 

of Tfl+l is three orders of magnitude less 
than the along-track position component 
of t n+ i. 


where Xq-i-i is the along-track state noise per 
unit time, D is a diagonal matrix with diagonal 
elements (1, .5, .2, 10“ 3 , .5(10)- 3 , .2(10)- 3 ), 
and ^(Tn+ 1) is a rotation matrix which trans- 
forms satellite state from an along-track, cross- 
track and radial coordinate set at time % to a 
standard inertial coordinate set. From Equation 
10 we have 


APn+1 “ (<Jx (Tn+1-Tn)^ (Tn+1 ) DDltyT (Tjh-x) (11) 

where 0 X is the standard deviation of . The 
value of 0 X can be chosen based on a priori con- 
siderations or it can be chosen adaptively from 
previous residuals. Adaptive algorithms based 
on a maximum likelihood principal [8], or on a 
minimum variance principal [9] are available and 
can be implemented on small computers. 

A more empirical approach to the problem is 
to account for the effect of errors in the prop- 
agation model in terms of an exponential process. 


P N+1,0 


= (<KT N+1 ,T N )P N <()T(T N+1 T N )) 
exp ^(Tfl+l- T N)/gN+l) 


( 12 ) 


The value of g^+i can then be chosen by adaptive 
techniques. 

Both the exponential approach of Equation 12 
and the additive approach implied by Equations 
11 and 4 were implemented for the real and simu- 
lated data reductions discussed in this paper. A 
totally automated adaptive procedure was not 
implemented. Instead, the data reductions were 
performed for numerous values of both 0 X of Equa- 
tion 11 and gjH-1 of Equation 12 and the values 
which maximized performance were chosen. 


Data Simulations 


In order to test the adequacy of the proposed 
procedure for data processing, the following 
simulations were performed. 

(i) Simulated data (see Figure 1) of range 
rate were generated for two stations (Portsmouth, 
N.H. and San Francisco, Ca.) relative to a satel- 
lite with the orbital characteristics of the 
TIROS-N satellite (see Table 1). The number of 
passes during which the satellite observed the 
transmitting stations is shown in Table 2. The 
satellite orbit was computed using a Goddard Earth 

Hodel gravity field [10] truncated to 13 x 13. 

The simulated measurements were corrupted by the 
addition of random gaussian noise with a stan- 
dard deviation of 0.5 meters /second. The 
choice of a 13 x 13 field was predicated on the 
assumption that the difference between the 
13 x 13 field and that of the real earth would 
be negligible in comparison with the difference 
between the 13 x 13 field and the field used in 
the subsequent orbit recovery. 
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Figure L> Data Simulation Configuration 
Two transmitting stations: 

(1) San Francisco, Ca 

(2) Portsmouth, NH 

TIROS-N IN NEAR POLAR ORBIT WITH A PERIOD 
OF 100 MIN. 


TABLE 1. Orbital Characteristics of TIROS-N 

Semi-Major Axis: 7200 KM 

Eccentricity: 0.000002 

Inclination: 98.7° 

Per iod : 100 Min . 


TABLE 2. Number of Observed' Passes In The 
Data Simulations 



0-8 Hr. 

8-16 Hr. 

16-24 Hr. 

Total 

DAY 1 

1 

4 

3 

8 

DAY 2 

1 

4 

3 

8 

DAY 3 

1 

4 

3 

8 

DAY 4 

2 

4 

2 

8 

DAY 5 

2 

4 

3 

9 

DAY 6 

1 

5 

3 

9 

DAY 7 

1 

4 

3 

8 

DAY 8 

2 

4 

2 

8 

DAY 9 

2 

4 

2 

8 

DAY 10 

2 

4 

3 

9 


(it) Two different fields were used in the 
orbit recovery: field I, a field containing only 

the J 2 harmonic coefficient and field II, a 4x4 
field. The non- zero coefficients of fields I and 
IX differed slightly from the corresponding coef- 
ficients of the 13 x 13 field. The filtering 
algorithm described in the preceding section was 
implemented to estimate orbits. The resulting 
orbits were then compared with the orbit used in 
the data generation. The results are summarized 
in Figures 2 through 4. 

Figure 2 shows the results for the J* field, 
with the covariance matrix of the state noise, as 


defined in Equation 5, arbitrarily set to zero. 

As can be seen, the error, which initially is 
very large, has been substantially reduced by the 
end of the first day. The error then starts 
growing. This is caused by force model errors 
which are not accounted for when the presence of 
state noise is ignored. The results of Figure 3 
were also obtained using a J 2 field but the state 
noise covariance matrix was obtained from Equation 
11 by setting 0 X = 4 M/HR. These results are a 
substantial improvement over those shown in Figure 
2. Figure 4 shows that more impressive results 
can be obtained with a 4 x 4 geopotential field 
with o x = 2 M/HR. It is also important to notice 
the long term stability of the filter when the 
effects of state noise are appropriately modeled. 
Comparable results were obtained using the ex- 
ponential model of Equation 12 with g n +l * 2 days. 
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Figure 2. 

Satellite Position Error 

Recovery Field Includes Only the J 2 Harmonic 
TIROS-N Simulated Data 

0 X ^= 0K/HR 



Satellite Position Error 

Recovery Field Includes Only the J 2 Harmonic 
TIROS-N Simulated Data 

0x - 4 M/HR 
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Figure 4. 

Satellite Position Error 
4x4 Recovery Field 
TIROS- N Simulated Data 
G x = 2 M/HR 
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Figure 5, 

Satellite Position Error 
4x4 Recovery Field 
NIMBUS-6 Doppler Data 
O x = 0.5 M/HR 


Reduction of NIMBUS-6 Doppler Data 

The orbit/ transmitter station' geometry for 
the NIMBUS-6 real data reduction experiment was 
similar to that used in the TIR0S-N simulations. 
The NIMBUS-6 orbit constants are given in Table 
3. Data from two transmitter stations were used 
in the NIMBUS- 6 orbit recovery, one station being 
located at Fairbanks, Alaska and the other one 
at Goddard Space Flight Center, Maryland. The 
nominal transmission frequency was 401.2 MHz 
for both stations. A gravity field complete 
through order 4 was used In the real data reduc- 
tion. Data for the first ten days of June, 1977 
were processed using the previously described 
recursive filter. 

The NIMBUS-6, orbit, determined as part of 
this study, was compared with an independently 
derived orbit based on Mini- track observations 
processed at GSFC, which has a claimed accuracy 
of about 500 meters. The orbit position dif- 
ferences shown in Figure 5 are less than 1 km 
RMS from the third day onward. As in the case 
of the simulated results, the filter’s perform- 
ance displays both accuracy and long term sta- 
bility. The number of passes of Doppler data 
used in the recovery is shown in Table 4. The 
noise in the Doppler data was approximately 1 m/s 
RMS. Data corresponding to elevation angles 
below 5 degrees were ignored. The data were not 
corrected for atmospheric refraction. The state 
noise covariance was chosen from Equation 11 
with a x * 0.5 M/HR. This constant was selected 
after a few trial riins. No attempt was made to 
fine tune the system. Comparable results were 
obtained with the exponential model with 
Sn+l = 5 days. 

In reducing the NIMBUS data, the oscillator 
bias and the bias drift rate were recovered for 
each transmitter* For the bias model constants 
actually used, the effective unmodeled drift rate 
contribution to the bias is shown in Table 5* 


TABLE 3. Orbital Characteristics of NIMBUS-6 
Semi-Major Axis: 7490 KM. 

Eccentricity: .001 

Inclination: 99.9° 

Period: 107 MIN. 


TABLE 4. Number of Passes Used in the Real 
Data Reduction 



0-6 HR 

6-12 HR 

12-18 HR 

18-24 HR 

Total 

DAY 1 

0 

2 

3 

1 

6 

DAY 2 

1 

1 

3 

1 

6 

DAY 3 

0 

2 

2 

1 

5 

DAY 4 

1 

3 

3 

1 

8 

DAY 5 

2 

1 

4 

3 

10 

DAY 6 

1 

4 

3 

2 

10 

DAY 7 

2 

1 

2 

2 

7 

DAY 8 

i 

2 

3 

3 

9 

DAY 9 

0 

3 

1 

1 

5 

DAY 10 

0 

0 

4 

1 

5 


TABLE 5. Growth of Bias Uncertainty as a 
Function of Time 

Time From Last 
Determination 
of Bias (Hours) 

Standard Deviation 
of Bias Uncertainty 
(M/S) 

1 

0.5 

2 

1.4 

4 

4.0 

8 

11.3 
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Conclusions 

For a class of satellite applications which 
involve a one-way Doppler positioning system, 
on-board orbit determination would be advanta- 
geous. This paper has presented a proposal for 
on-board orbit determination which would rely 
entirely on reference Doppler beacons located 
within the Continental United States. 

To demonstrate the feasibility of this 
concept it is necessary to develop algorithms 
which are compatible with small computing 
machines and sufficiently accurate to maintain 
an adequate satellite ephemeris . The simula- 
tions and real data reductions discussed in this 
paper suggest that an on-board computer equipped 
with a recursive filter with a fading memory 
to account for dynamic modeling errors, and two 
reference beacons, are adequate for the task. 

The filter should include a gravity field 
model with geopotential coefficients complete 
to degree and order 4, The fading memory 
characteristics of the filter can be set in 
advance based on a priori simulations and real 
data reductions. An alternative possibility is 
to equip the filter with a feedback mechanism by 
which the memory would become a function of 
previous residuals. Several such procedures 
have been discussed in the literature. The 
core storage requirements of the resulting 
algorithms should be compatible with on-board 
computers. 
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